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Abstract 

The TTT-TTTT spectrum and various angular distributions in T(3S') — > T(IS') vrvr 
are studied including the effects of the vrvr phase shift in the I = L = chan- 
nel using the lowest order amplitude in the pion momentum expansion. Our 
results are compared with the recent CLEO data, and we find good agreement 
except for the cos0* distributions. We argue that the cos 6* distribution, con- 
trary to other distributions, is sensitive to the higher order corrections in the 
pion momentum expansion. This argument is supported by using an ansatz 
for the amplitude which is of higher order in the pion momentum expansion 
and still satisfies the soft pion theorem. 

PACS Numbers : 14.40. Gx, 13.25.+m, 13.40.Hq 
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Typeset using REVT^ 



I. INTRODUCTION 



The peculiar double peaks in the m^^r spectrum in T{3S) T(IS') rni have been 
lacking proper understanding. Although several suggestions have been made, most of them 
reproduce the tttt spectrum only with limited success ^ . In the two recent works by us, we 
approached this problem in two different ways. In the first work we assumed that the 
most general form of the amplitude for T(35') — > T(IS') tttt is 

M = Ao + B0E1E2 + Coml) e ■ e + Do (p ■ ep ■ e + p ■ e p ■ e)] , (1) 

in the lowest order in the pion momentum expansion. Here, p,p' are the four-momenta of 
the final pions, Ei and E2 are their energies, and e and e are the polarization vectors of 
the initial and the final T's, respectively. In Ref. 0, we found three sets of parameters (PO, 
PI and P2) by minimizing the x^- Fo^' these three sets of parameters, we predicted various 
angular distributions which should be checked against experiments. 

In the second work [^], we assumed that (i) the QCD multipole expansion is applicable 
to T{3S) —>■ T(15') TTTT and (ii) T(35') has an admixture of a D— wave component : 

\T{3S)) = cos (l)\3S) + sin (f)\D). (2) 

The amplitude for T{3S) T(IS') nn then depends on three independent parameters, 
and we found two fits which correspond to PI and P2 of Ref. 0, respectively. We then 
explored the consequences of our assumptions on other hadronic and radiative transitions of 
T{3S) into the lower level bottomonia. In Ref. |Q, we again assumed that the lowest order 
expansion in the pion momenta 

^ {nnr^lO) = {q' + ml) (3) 

be valid through the whole range of m^rTr = 

Although these two approaches fit the 171^^^^ spectrum, they still leave room for theoretical 
improvement in two aspects. First of all, it is well known from the analysis of the vrvr phase 
shift that the dipion system in / = L = experiences strong final state interactions [^. 



Since the dipion system in T{3S) T{1S) vrvr are in either J = L = OorJ = 0, L = 2 
state, one should properly take into account of the tttt phase shift due to the final state 
interactions in the I = L = dipion system. 

Secondly, the validity of (|lj) or (|) in T{3S) T(15') vrvr is rather unclear, since the 
available m^rTr in T{3S) T{1S) tttt is large, 2771^^ < m^^ < (m^ — ruf) = 895 MeV. The 
amplitude ([l|) with Bq = Dq = gives a good description for the m^^r spectrum in T{2S) — >• 
T(IS') TTTT, where m^rTr < 563 MeV. For higher value of m,r7r, we simply assume that it is valid 
and explore its consequences on various spectra of decay products in T{3S) T(IS') tttt. If 
any of the predicted angular distributions based on the amplitude (|I]) including the tttt phase 
shift does not agree with the experimental measurements, it would signal the importance of 
higher order terms in the pion momenta which have been neglected in (0). 

In the present work, we follow the approach of Ref. and include the phase shift of the 
TTTT system in the I = L = channel ((5o(q'^)) using the data available in the literature 0. 
The phase shift for the 1 = D—wave tttt system (52(9^)) is tiny enough to be neglected for 
the whole range of m^^ [0. 

In Sec. p.1 A| , we decompose the amplitude (p into the tttt 5*— and D— waves, and incor- 
porate the phase shift Solq"^). It is found that the importance of the phase shift due to the 
final state interactions is most prominent in the cos^* distributions, but the effect is only 
moderate. Then, our results are compared with the recent data from CLEO in Sec. |IIB| . In 
Sec. m, we generalize the amplitude ([l|) (in case of D = 0) to include higher order terms in 
g^, and discuss some general aspects of the various angular distributions. It is also pointed 
out that one can extract the difference between tttt phase shifts of the S— and D— waves 
in the 1 = channel, — ^2, by measuring the joint distribution d^r/dm,r7rdcos6'*. The 
results are summarized in Sec. fTV] . 



II. EFFECTS OF THE vrvr PHASE SHIFT 
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A. Predictions on the cos 9-,^* distributions 



In this work, we use a modified form of tfie amplitude (|T]). First of all, we ignore the recoil 
of the final T(IS') in the amplitude, and make the D— term proportinal to the symmetric 
traceless part : 



M = A 



+ BE1E2 + Cm^l e ■ e' + |p ■ ep ■ e + p ■ ep ■ t - '^p- pe ■ e'j 



(4) 



One can find relations between parameters in (|l|) and using 

p-f^ = E,E2-^ is-2ml). 

It should be emphasized that our amphtudes (P and (0) satisfy Adler's condition by con- 
struction. 

Now, we decompose the above amplitude into the vrvr S"— wave and D— wave in order to 
take into account the phase shift of the vrvr system in the I = L = state. The S— and 
D— waves have the following tensor structures 0] : 



V,^= fr,{q') [cos' ei--] 5,,+go{q' 



riTj - -{qtqj + q^S^j)f3f 



(5) 



where r^= p^ — p^. The amplitude (H) can be decomposed into the fs,D,gs,D form factors 
as follows : 



fs = q' + Cml + 



B 



{E, + E,f - - <f^f 



- Dq'--D 



-2 a*2 / 2 I -2 



= i £> (1 - i/j;') 



90 = -5 D. 

Multiplying the 5*— wave amplitude by the phase shift 5o(g^), we get 



(6) 
(7) 
(8) 
(9) 



M 



(10) 



In case the final T(15') is not reconstructed, summations over polarizations of the initial 
and final T's are done with 

^iiij = {5ij - ZiZj) , (11) 

E^i^; = ^^.- (12) 

The fact that the initial T(35') is transversely polarized with respect to the beam directions 
(taken along the z— direction) has been taken into account in (|ll|). This fact is very useful 
to test the existence of the D— term by measuring the polar angle distributions of the final 
T(IS') (or, equivalently, the nir system as a whole) and/or of a muon emerging from the 
muonic decay of the final T(15'). If one tags the muonic decay of the final T(15'), the 
polarization sum over the final T(15') ([121) should be replaced by 

Ee'ii'j = {6ij -Uj), (13) 

where / is the three dimensional unit vector along the direction of a muon in the rest frame 
of the initial T(35'). For D = 0, one gets 

dN /d cos 6i ~ (1 + cos^ 6i), 

where cos6'; = I ■ z. 

Using amplitude (|10|), we fit the m^rTr spectrum by minimizing x^- The best fit is given 
by three sets of solutions, PO, PI and P2 (see Table 1), which are essentially the same 
as the ones given in Ref. 0. Other angular distributions can be obtained by numerical 
integrations as in 0. In Fig. |I| (a) and (b), we show the cos^* distributions of vr^ in the 
rest frame of the dipion system for PO and PI (P2), where 6* = 0° is along the direction of 
the dipion system as a whole in the rest frame of initial T(35'). For comparison, we show 
the corresponding plots with the phase shift neglected {i.e. 5o(g^) = 0) in Fig. 2 (a) and (b). 
The phase shift moderately changes the cos 6'* distributions, but the overall effects may be 
hardly discernible in the experiment. 

Before continuing on to the next section, we discuss how much the results obtained in 
Ref. will change when we incorporate the vrvr phase shift and possible corrections to (3) in 
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higher orders in the pion momentum expansion. In this case, it suffices to resort to the QCD 
multipole expansion by our assumptions. The relevant matrix element for (7r7r|Gj^j^G'"^'^|0) 
has been obtained by Donoghue et al. using the dispersive approach in conjunction with 
the chiral symmetry relations imposed in the chiral limit 0. The result is that the above 
matrix element, (7r7r|G'^j^G"^^'^|0), is dominated by (7r7r|^^|0), and that (3) remains essentially 
unchanged up to m.„T, ~ 0.9 GeV, once it is regarded as the modulus of 02-k- The phase shift 
is given by 5Q{q^). In short, one only has to write (3) as 

(7r7r|^;:|0) = (g2 + m^)e'^«(<^'). (14) 

Since the phase shift ^Q^q^) does not affect the m^r^ spectrum in T(35') T(IS') vrvr, our 
results derived in Ref. do not change at all. In particular, the predictions for T(3S') — >■ 
T(IS') + 7] remain the same, which excludes the fit PI ||^. 

B. Comparisons with the data 

Recently, the CLEO collaboration released a new set of data on hadronic transitions in 
T(3S') decays 0. Their results on T(3S') T(IS') tttt can be summarized as follows : 

(i) the cos Of distribution is flat for the whole range of cos6f. 

(ii) the cos^TT* distributions can be fitted by (a + 6cos6'^*), with a = (1.24 ± 0.06) 
and b= (-0.49 ± 0.13). 

(iii) the cos6'; distribution is consistent with (1 + cos 9f) for < | cos^^;| < 0.7. 

Let us discuss the implication of each statement above to our fits, P0-P2. Statement (i) 
excludes both PI and P2, since these two lead to quadratic functions of cosOj. Statement 
(iii) also partly supports this conclusion, since (iii) implies D = Oin the amplitude (1). Thus, 
(i) and (iii) select PO as the final candidate. However, the cos 6*^* distribution shown in Fig. |I| 
(a) does not agree with statement (ii) from CLEO. We do not interpret this disagreement 
of the CLEO data with our prediction on the cos ^* distribution as a general failure of our 
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approach based on the matrix element satisfying the soft pion theorem. We rather regard 
it as an indication that amphtude (1) needs to be modified to include higher order terms in 
the pion momentum expansion. This will be illustrated in the next section with a modified 
amplitude for T(35') T(IS') tttt. 



III. MORE ON THE AMPLITUDE WITH D = 



In this section, we consider the case D = in more detail, including possible higher order 
terms in in the 5— wave tttt amplitude, fs- It will be shown that the cos^* distribution 
is sensitive to such higher order terms in contrary to other distributions. 

Let us write the amplitude for T{3S) T(IS') tttt as 



M{T{3S) T{1S) tttt) = A 



e ■ e 



(15) 



where fs and fc satisfy the soft pion theorem. The explicit forms of fs^n for the lowest 
order amplitude @ can be read off from with D = 0. Therefore the differential cross 

section for e^e^ T{3S) T(IS') tttt mifi~^fi~ is given by 



di^T oc dm-^T^d cos 9l.d cos 9i m,r7r|i?| PI 1 + cos^^/ 

fs + fl (cos^ - ^) ' + VsfD COS 5 {cos' e:-^- 



(16) 



where /?* is the velocity of a pion in the tttt rest frame and 9i is the angle between a muon 
and the e^e~ beam in the rest frame of T(3S'). 

Integrating the partial distribution ([16|) over appropriate variables, one gets 

dV 



dcos Of 
dV 



oc 1, (fiat distribution), 
oc (1 + cos^6'i). 



(17) 
(18) 



d cos 9i 

These two distributions are independent of the dependence of the form factors, /^(g^) 
and foiff)-, as well as of the vrvr phase shift. And, these results are consistent with the recent 
report from CLEO. 



On the other hand, the cos6'* distribution is sensitive to the actual forms of fs{<f) and 
foiff)-, and to the vrvr phase shift. In principle, there are many possible terms to the next 
order in the pion momentum expansion. Instead of writing down all possible terms and 
fitting the m^jr spectrum as in Ref. H, we take the following amplitude for illustration : 



l + C 



EI + E2 



B 
+ 4 



{E, + E,f--\pf\'(3f 



(19) 



with the same fniQ'^) as before. This amplitude has three parameters, A,B and C, and 
satisfies the soft pion theorem like (|^). 

By fit to the 171^^^^ spectrum, we found another fit (we will call it P3) with x^/d.o.f. = 
11.0/7 (see Fig. ^ (a)). The corresponding values of A, B, C are given in the last column of 
Table 1. This amplitude predicts the distributions, (|T7p and (|TB|). The cos6'* distribution 
for P3 shown in Fig. ^ (b) differs a lot from that for PO in Fig. |l] (a), and gets much 
closer to the observed data. The lesson from this example is that once we adopt (p!5[) as the 
amplitude for T(3S') T(IS') tttt, we predict (i) the flat cos6'/ distribution, (ii) (H-cos^6'i) 
distribution for the polar angle of a muon, independent of actual forms of fs, fn and 5o(g^). 
This is not the case for cos 6** distribution and thus cannot be reliably calculated unless 
they are known. The actual functional forms of fs, fn and 5o(q'^) can be extracted from the 
measurement of the joint distribution, (f T / drnT^-^d cos 9* , as one can derive from (|T6|) : 



dm.^Tjd cos 6^ 
oc m^^\q\Pl 

= "Co(g2) + C2(g2) cos^ + C,{q^) cos^ 61 



fs + foi cos^ e:-^j + 2fsfD COS 60 [cos' 91 g 



(20) 
(21) 



For each 771.,^^^ bin, one can measure the cos^* distribution. This determines Cj(g^)'s, and 
in turn, three unknowns, fs.fo and 5o(q'^)- In particular, the decay T(3S') T(IS') vrvr 
can be a source of the 5— wave vrvr phase shift for the whole elastic region, 2m^ < m^rTr < 
(m, — m/) = 895 MeV |T0|. This may be important, since the existing data on the tttt phase 



shift between itlk < ^t^tttt < 600 MeV are rather poor in statistics and one has to make some 
extrapolation |Tl|. 
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IV. CONCLUSION 



Concluding, we reanalyzed the T(3S') T{1S) tttt decay using the most general matrix 
element in the lowest order in the pion momentum expansion, including the final state 
interactions of the vrvr system in the I = L = channel. The vrvr phase shift changes the 
cos 6'* distributions moderately. (Compare Figs. 1 (a), (b) with Figs. 2 (a), (b).) Compared 
with the recent data from CLEO, PO is selected, but the cos6'* distribution does not agree. 
In Sec. |TT|, we argued that this distribution is sensitive to possible higher order corrections in 
the pion momentum expansion. As an illustration, we used a new ansatz for the vrvr S*— wave 
amplitude, ([19|) , which is of higher order in the pion momentum expansion, and satisfies 
Adler's condition. This amplitude could fit the m^rTr spectrum (Fig. 3 (a)). The resulting 
cos 6'* distribution (Fig. 3 (b)) is different from Fig, 1 (a), and gets closer qualitatively to 
the measured distribution although not quantitatively. It would be more proper to do the 
partial wave analysis with the S— and D— waves, and find out the form factors {fs, fo) 
and the phase shift ((5o(g^)) in ([T5|) from the joint distributions in m^^r and cos 6** using 
(20) and (21). Since our explanation involves both S— and wave tttt systems, we may 
be able to find out the S'— wave phase shift (or 8^ — 82, more precisely) from the decay 
T(3S') T(IS') TTTT by measuring various joint distributions. Since the available m^^r is 
below the KK threshold, this decay may provide information on the phase shift over the 
whole elastic region, especially for mx < "n^vrTr < 600 MeV, where the current data are rather 
poor in statistics. 
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FIGURES 

FIG. 1. The cos^* distributions of 7r+ in the rest frame of the dipion system with the phase 
shift included : (a) for PO and (b) for PI (and P2). 6* = 0° corresponds to the direction of the 
dipion system as a whole in the rest frame of the initial T(3S'). The dotted and the dashed curves 
correspond to the low and the high mTTTr regions, and the solid one represents the sum of the two. 

FIG. 2. The cos^* distributions of 7r+ in the rest frame of the dipion system with the phase 
shift neglected : (a) for PO and (b) for PI (and P2). 9* = 0° corresponds to the direction of the 
dipion system as a whole in the rest frame of the initial T(3S'). The dotted and the dashed curves 
correspond to the low and the high mTTTr regions, and the solid one represents the sum of the two. 

FIG. 3. (a) The mTTTr spectrum and (b) the cos9* distributions for the amplitude, (15), where 
fs is given by (19) and fn is given by (8) with D = 0. For (b), the dotted and the dashed curves 
correspond to the low and the high 171^^^^ regions, and the solid one represents the sum of the two, 
respectively. 
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TABLES 

TABLE I. Three sets of parameters for the amphtude (10) and one set for the ampUtude (19) 

giving the best fit to the m7r,r spectrum. PO and P3 correspond to the constrained fit with 
D = 0. The parameter A is the overall normalization. 



Fits PO PI ( P2 ) P3 

A 10.5 ±0.6 6.7 ±2.1 163.6 ± 55.3 

B -6.4 ±0.7 -3.7 ±3.9 -0.140 ±0.035 

C 36.2 ±6.2 7.8 ±37.2 -0.154 ±0.001 

D 0.0 (fixed) Tl.7±0.7 0.0 (fixed) 

X^/d.o.f. 12.0 / 7 11.1 / 6 11.0 / 7 

C.L. 10.2 % 8.4 % 14.0 % 
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